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STRUCTURED  ELECTRODE  INTERFACES 


LESSER  BLUM 

Department  of  Physics,  PO  Box  At,  University  of 
Pueno  Rico,  Rio  Piedras,  Puerto  Rico,  00931 

Abstract 

A  general  discussion  of  the  structure  of  the  charged  interface  when 
the  surface  is  not  smooth  is  presented.  In  particular,  the  effect  of  a 
structured  electrode  is  presented. 

Introduction 

- The  electrochemical  interface  offers  a  particularly  big  challenge 

to  the  scientist:  because  it  is  such  an  important  system  a  lot  has 
been  worked  on  it,  and  yet  because  of  the  intrinsic  technical 
difficulties  associated  with  non-extra  high-vacuum  surface  physics, 
very  little  direct  structural  evidence  is  available  for  the  active 
electrode. 

There  is  an  increasing  body  of  evidence  coming  from  optical 

spectroscopy-Techniques  such  as  Surface  Enhanced  Raman  effect 

_S  c 

(SERS),pT  Surface  Infrared  spectroscopy  (SIRS)  ,  and  more 

recently,  second  harmonic  generation  (SHG)  are  fine  examples  of 

/ 

diagnostics  for  structural  changes  (changes  in  bonding  arrangements) 

but  do  not  provide  unequivocal  structural  information:  The  shifts  and 
intensities  of  the  spectral  lives  depend  on  local  field  as  well 

as  local  structure.  ,  <  -/ 

Quite  interesting  information  can  be  obtained  from  ex-situ 
experiments,  in  which  the  electrodes  are  examined  by  surface  sensitive 
techniques  [4],  (LEED,  for  example  in  extra  high  vacuum.  As  has  happened  in 
the  science  of  solid  surfaces,  we  learn  that  the  structure  of 
electrodes  is  rich  and  complicated.  For  example,  there  may  be  ordered 
2-dimensional  phases  formed  by  species  present  at  the  interface.  Take, 
for  example  the  case  of  Pb  deposited  on  silver  (111).  This  is  a  case 
of  under  potential  deposition.  A  mono  layer  of  Pb  is  deposited 
much  before  bulk  lead  is.  Kolb,[5]  and  co  workers  have  performed 


ex-situ  studies  of  this  system  and  shown  that  there  can  be  two  phases 
of  Pb.  One  is  a  dense,  hexagonal  phase  which,  however,  is 
incommensurate  with  the  underlying  Ag(lll)  surface.  A  low  density 
phase  (1/3  of  full  coverage)  is  also  observed  in  these  experiments. 

At  intermediate  coverings  a  mixture  of  the  two  phases  is  observed. 


The  question  is:  what  happens  in  the  active  electrode?  It  is 
quite  fortunate  that,  direct,  x-Ray  observations  can  (and  recently 
have[6]  been  made  on  this  system. 

The  results  indicate  that: 

a)  at  high  voltages  (-.53  to  -  1  v)  only  the  dense  phase  exists. 

b)  at  lower  voltages,  (-.38,  -.39  v)  where  the  shoulder  in  the 
voltamogram  is,  the  low  density ,commensurate  -J3  x  J3  R  30°  phase  is 
not  observed. 

c)  furthermore,  at  intermediate  coverage,  only  the  dense  phase  exists 
in  the  form  of  clusters.  These  clusters  grow  in  size,  (and  number) 
shrink  and  twist  slightly  during  the  deposition,  which  takes  place 
between  -.375  and  -.525  v.  from  x-Ray  diffraction  experiments[7])  .The 
film  seems  to  be  stable  between  -.53  V  and  -1  .V. 

The  question  now  is:  how  do  we  formulate  a  statistical  mechanical 
theory  of  the  double  layer  (in  this  case,  the  inner  Helmholtz  layer)  in 
which  the  structure  of  the  solid  surfaced  is  included  -  clearly  this  is 
an  old  problem  that  has  been  discussed  by  many  distinguished 
researchers  like  Frumkin,  Levine,  Fawcett,  just  to  mention  a  few  names. 
But  their  approach  was  necessarily  more  empirical  than  fundamental.  Our 
point  of  view  is  going  to  the  different.  We  will  propose  simple, 
schematic  models  of  the  structured  interfaces,  about  which  we  can  make 


precise  statements.  (It  is  similar  to  the  drunkard  who  lost  the  key  to 

his  home  somewhere,  but  looks  for  it  under  a  light  pole,  because 

it  is  the  only  place  where  he  can  see).  So,  let  us  consider  a  smooth 


electrode  again  (one  sided,  meaning  a  hard  impenetrable  surface,  or  two _ 

sided,  which  consists  of  a  surface  separating  two  conducting  media), 

and  let  us  decorate  it  with  sticky  sites[8]  .  nTTr  tab 


exp  [-|3u(r)]=  1  +  X5(r-a') 


«TIS  GRAtl 
OTIC  TAB 
Unannounued 
(1  )Ju3tif laatloi 


where  8  (r)  is  Dirac’s  delta  function,  (3  is  Boltzmamn’s  thermal  - 

factor,  and  X  is  the  sticking  probability.  Again  we  start  at  the  light 
pole,  like  the  drunkard:  Our  model  may  not  be  very  close  to  the  real  at 


Speoial 


□  □ 


interface,  but  it  has  been  simplified  so  that  we  may  be  able  to 
understand  it.  Let  us  now  turn  to  an  even  simpler  version 

(Jancovici[9]):  In  two  dimensions,  the  one  component  plasma  (OCP),  that 
is,  the  system  in  which  the  positive  ions,  for  example  have  been  smeared 
to  become  a  uniform  neutralizing  background,  and  the  negative  ions  are 
point  charges.  The  negative  ions  interact  through  the  potential. 

u(r..)  =  -e^  In  tr  .l 

IJ  *  IJ  ' 

For  a  specific  temperature  T,  for  which 

^W/kT=2 

The  statistical  mechanics  of  the  OCP  is  completely  solvable  .  This 
means  that  we  can  use  this  model  as  a  tool  (of  "experimental 

mathematics")  to  test  exact  sum  rules,  valid  for  the  more  realistic  3 

dimensional  models.  Very  briefly,  to  see  why  this  model  is 

exactlysolvable,  consider  the  configurational  canonical  partition 

function  of  a  system  of  N  negative  charges 


(2  ) 

(  3  ) 


Z,(X)  =- 


NT 


dr^..dr^  exp  (-P([  v.)l 


1 

NT 


drj..dr^  exp  v  J]  exp  [  (3([  e^  In  |r,  J)] 


1 

NT 


dr^.Ar^  exp  [-P([  v.)]  |n(r-  r. 


(4) 


Here,  v,  =  v.(r.)  is  the  interaction  of  the  charge  with  the  background. 
In  the  simplest  case,  it  will  be  the  uniform  charge  density,  in  which 
case  v.oe  |rj^ 

There  seems  to  be  another  two  ingredients  that  are  necessary  to 
obtain  explicit  solutions  for  these  models!  10] 

a)  that  the  product  in  equation  (4)  can  be  written  as  a  determinant  of 
dimension  NxN. 

b)  that  the  determinant  can  be  diagonalized  for  arbitrary  values  of  N. 
This  diagonalization  has  been  achieved  for  backgrounds  which  mimic  one 
and  two  sided  electrodes,  and,  in  general  for  arbitrary  forms  of 

V 


v,(x)  (Here  r=  x,y). 

More  recently  it  has  been  shown[ll]  that  the  case  of  a  structured 
interface  of  a  two  dimensional  OCP  can  be  soKv''  explicitly.  In  the 
simplest  case,  v=  x+  V  and  V  is  defined 

I  i  ,st  SI 

M 

exp  [-PV^^(r)]=  1  +  A.  ^  8(y-yJ&(\) 

m=  1 

(5  ) 

where  y  is  the  position  of  the  m^^  sticky  site.  All  the  sticky  sites 

Tn 

are  located  at  the  straight  line  defined  by  x=0,  and  the  system  is 
otherwise  uniform.  M  is  the  total  number  of  sites  for  a  finite  length  L 
of  the  line. 

0)=M/L  (6  ) 


is  the  density  of  the  sites  on  the  line.  The  excess  firee  energy  Af* =AF/L 
due  to  the  presence  of  the  line  of  adsorbing  sites  is  then 


AF=-  kT  In 


(7  ) 


Af‘=AF/L=  -kT  0)  JdC  ln(l+Xp  0^(C,t)] 


where 


t=7cp/2co^ 


p  is  the  density  of  charges  in  the  bulk  system,  and  is  the 

Jacobi  theta  function  defined  by[12] 


0  (c.t)=  y 

^  m=-« 


g-tm^+2Jtim^ 


(9) 


We  can  also  compute  other  properties  of  this  model.  The  fraction  of 
occupied  sites  is  (the  adsorption  isotherm!). 


n 

« 


03(C.t) 

[1+xp  ejCC.t)] 


and  the  density  profile  p(r)  is,  using 

M 

P(r)=y(r)[l+  X  ^  S(r-rJ] 

m=  I 


with 


y(r)=  p  -  p^  e-^Pl^l' 


F(C,z,t) 

[1+Xp  e3(C.t)] 


2F(^,z,t)  =  e3(C+t£Oz/7C,t)  03(^+tO)Z*/jt,t) 

+  ta)z/)t,t)  S^CC-tcoz  /7t,t) 


and  z  =  X  +  iy;  z*  =  X  -  iy. 


Consider  now  the  case  of  the  one  sided  interface  (a)  or 
interface  (b).Both  of  these  cases  were  recently  discussed  by 
We  follow  her  analysis  closely.  For  case  (a):In  that  case  we  get 

i+o/co 

Af=  -kT  (0  [  dC  ln[l+Xp  Hj(0] 
o/co 


n 

a 


i+o/co 


=  J  dC 

a/o) 


H  ,  (O 

[1+Xp  Hj(C)] 


where 


H,(C)  =  S  (C+n) 


Sj(u)=  2Vn/Vt 


JC^uVt 


(10) 

(11) 

(12) 

(13) 

(14) 

two  sided 
Cornu  [13]. 

(15) 

(16) 

(17) 

(18) 


l+OfTtu/t) 


du 


<I>(u)=  INk  fdv  e*'^ 

J  n 


(19) 


and  the  single  particle  density  is  (regular  part) 


l+O/O) 


y(r)=  p  e  ^'^P*  I  I  S  ^  ( C+n ) 

o/co  " 


(20) 


,  27tO)z  (n-t-^) 


Xp 


[lUp  Hj(0] 


[  Sj(C+m) 


J 


For  the  double-sided  model  of  the  interface  (b)  we  have 
Excess  free  energy 

1 

Af=  -kT  0)  J  dC  In  [  1+Xp  H^(C)  j  (21) 

0 


Adsorption  isotherm  (fraction  of  occupied  sites) 


n  = 

a 


I  dc 


Xp^H2(C) 


(22) 


Single  particle  density 


y(r)=p.e-^>'P.'%-2(^-^)  J  I  S,(C^n)  ' 


(23) 


g27CO)z*(n+C)  .  ^  j-  ^  g27CO)z’(m+; ) 


l+Xp^H^CO 


where 


Sfi+n) 


(24) 


-K^U^/t 


S^(\i)=  2Vk/A 


2  2 

erfc(-7cuv't)  +m  e*^  ^  erfc(mreuv't) 


(25) 


with 


erfc(u)=l-<I>(u) 


(26) 


m= 


(27) 


Nest  question  is:  can  we  draw  conclusions  from  the  two  dimensional 
models  that  will  be  also  valid  for  the  three  dimension  case? 

Consider  now  the  case  of  an  impermeable  wall  in  3d,  decorated  with 
an  array  of  sticky  sites  -  we  use  the  definition 

Pi(n=yj(l)  exp[-pv,  (r)]= 

=y,(l)[  1  +  5(z)  6(R-njaj  ^283)  ]  ™ 

"."2 

where  p.(l)  is  the  density  of  species  i,  r=  (x,y,z),  R  =  (y,z),  and  the 
sum  n^.n^  is  over  all  the  sites  on  the  surface.  Then,  we  can  derive  a 
contact  theorem  that  will  relate  the  contact  density  to  the  charge  on 
the  wall  and  thermodynamic  properties  of  the  fluid  [14]  .  Assume  that 
the  sticky  sites  form  a  regular  lattice,  and  that  the  area  of  the  unit 


cell  in  that  lattice  is  A.  From  (28),  we  write 


^  ^  -Pv(l) 

-kT  Vjp,(l)=  p.(l)  VjUXD  -  yil)  Vj  e  ‘ 


^^[d2p,  X12)  VjU_.(12) 


where  v,(l)  is  the  sticky  potential  (1)  alone  and  all  the  other 
contributions  to  the  single  ion  potential  are  included  in  U.(Xj)  and 
u.  (12)  is  the  two  body  potential.  This  equation  is  very  complicated 
because  it  is  three  dimensional.  If  we  now  take  averages  over  the  unit 
cell,  that  is,  we  define 


y  (0)  =  1/A 


dx^dy^  p.(Xj,yj, 
*^cell 


then  we  get  by  integration 


I  1  ’ 

2 

where  is  the  bulk  pressure,  e  is  the  dielectric  constant,  <E^>q 
is  the  mean  square  bare  (vacuum)  field  at  the  undecorated  interface. 

This  contact  theorem  is  verified  by  the  exactly  solvable  model 

Let’s  then  try  to  give  a  more  quantitative  treatment  of  this 
model[15]  .  We  start  by  a  model  in  which  we  have  a  hard  sphere  fluid  in 
contact  with  a  wall  decorated  with  sticky  sites.  Let  us  call  the 
hamiltonian  of  the  N  particle  system  without  the  sticky  sites,  and 
the  sticky  potential 


H=T  v(l) 

’  .^1  ‘ 


Then 


Z(0)=-j5rT  Jdr'^  e  ® 


1  r..N 


Z(X)=-j5ry  Idr""  e 


(33) 


and  because  of  (1) 


Z(X)=- 


N! 


dr^  e 


-pH  N 


°  n  1+  8(zp5(R-R.) 


i  =  l 


(34) 


Divide  and  multiply  by  Z(0),  and  recall  the  definition  of  the  s-body 
distribution  functions,  then 


■ZHJT  -  L  —^r 


(35) 


with 


p°(Rj,...Rp  Z(0)N!  ® 


(36) 


=  g°(R,,...Rpnp!(R) 


where  p.(r)  is  the  singlet  distribution  function  for  the  undecorated 
model.  Using  standard  thermodynamics,  and  the  fact  that  X.  can  be 
identiHed  up  to  a  multiplicative  constant  (or,  equivalently  up  to  a 
shift  in  the  reference  chemical  potential)  with  the  fugacity  of  the 
adsorbed  particles,  we  get 

^  >"  [^] 

Af=AF/A 

S 

and  the  Gibbs  equation 

Ar=  r(X.)-r(0) 

where  r(X)  is  the  total  surface  excess  (decorated)  and  r(0)  the 
undecorated  one. 

Clearly  starting  with  (35),  we  have  formal  expressions  to  compute 
the  properties  of  the  interface  .  H  o  wever  we  need  to  know  the 
distribution  functions  of  the  undecorated  system  to  all  orders. 

We  did,  indeed,  know  this  function  for  the  2-D,  solvable  model.  But  we 
do  not  know  it  in  general.  So  we  have  to  make  approximations,  the 
simplest  of  which  is  Kirkwood’s  superposition  approximation  which  is 


(37) 

(38) 
(39) 

(40) 


If  we  now  introduce  the  potential  of  mean  force  between  panicles  in 
sites,  i,j  of  coordinates  r  =  (0,  R),  r=  (0,  R) 

I  1  j  j 


0),  =  -kT  In  g--{r  r) 

ij  1 } 


(42) 


Then  we  can  rewrite  (35)  as 


Z(X)  _ 

ZM 


All  5  ites 


s! 


cxpj^-Pj^  W2(r,r)j  (43) 


_  -  lies 

on  the  surface 


We  observe  now  that  this  is  just  isomorphic  to  the  grand  canonical 
partition  function  of  a  two  dimensional  lattice  gas  in  an  external 
field.  The  site-site  interaction  is  given  by  the  potential  of  mean 
force  (0  .  ,  and  the  external  field  by  a  potential 


V  =  -kT  In  p  (0) 

Using  the  Percus-Yevick  expressions  for  g^ftJ  )  and  for  p''(0) 

g"(o  )  =  (1  +ti/2)/(1  -  Ti)' 

p°(0)=  (1  +  2ti)  /(I  -  Ti)^ 

il=7t/6pa^ 


(44) 


(45) 


we  get,  using  mean  field  theory  (Bragg- Williams  approximation). 


n  = 

a 


n  c 

X  pj(0)|^g°(a)j 

n  c, 

1  +  ^  P”(0)[g^(a)j  I 


(47) 


where  c  is  the  number  of  nearest  neighbors  in  the  lattice.  This  is  an 
equation  for  the  fraction  of  occupied  sites  n^.  At  high  densities  and 
low  coupling  constants  there  is  a  phase  transition.  The  reason  for 
this  is  that  the  contact  probability  g2(<J)  inc’^eases  with  density, 
meaning  that  the  effective  pair  attraction  of  particles  in  occupied 
sites  also  increases  with  density.  The  effect  of  the  sticky  parameter  X 


(actually,  of  X  p°(0)  )  is  that  of  the  external  field,  which  has  the 

effect  of  attracting  indiscriminately  to  any  site  of  the  interface,  and 
therefore  has  a  disordering  effect. 

Very  recently  Caillol,  Levesque  and  Weis[16],  have  carried  out  a 
detailed  Monte  Carlo  simulation  of  a  system  that  is  similar  to  the 
sticky  surface  system.  The  major  difference  is  that  the  sticky  sites  are 
not  points,  but  rather  small  muffin  tins  (camemben  boxes)  of 

cylindrical  shape.  Any  ways  the  simulation  does  not  see  a  sharp 
transition,  as  predicted  by  the  theory.  However,  the  qualitative  effect 

of  clustering  is  clearly  noticeable. 

Furthermore,  the  size  of  the  clusters  increases  with  the  coverage 
n^,  and  furthermore,  they  are  not  really  commensurate  with  the  lattice 
of  adsorbing  sites.  This  is  due  to  the  fact  that  because  of  thermal 

rattling,  the  ’effective’  diameter  of  the  adsorbate  is  bigger 
than  o.  Therefore  the  clusters  are  not  commensurate  with  the  surface.  But 
certainly,  this  model  does  imitate  the  features  of  our  Pb/Ag(lll) 

interface,  at  least  qualitatively  in  the  sense  that  the  effect  on  the 
higher  coverage  is  to  decrease  the  distance  between  the  adsorbed 
particles,  while  the  size  of  the  clusters  increases.  To  see  a  sharp 
transition  we  probably  have  to  make  the  camembert  boxes  smaller,  or  even 

introduce  an  attractive  potential  between  the  adsorbed  particles.  This 

does  not  seem  to  be  the  case  in  the  computer  simulation,  or  the 
experimental  data  on  the  Pb/Ag(lll)  system. 
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